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ABSTRACT 

We present spectroscopic follow-up observations of galaxies from the Munich Near-Infrared 
Cluster Survey (MUNICS). MUNICS is a wide-field medium-deep A"'-band selected survey 
covering 1 square degree in the near-infrared K' and J pass-bands, and 0.35 square degrees 
in I, R, V, and - recently completed - B. The spectroscopic sample comprises observations 
of objects down to a limit of K' < 17.5 in five survey fields (0. 17 square degrees in total), and 
a sparsely selected deeper sample (if < 19.0) constructed in one of the survey patches (0.03 
square degrees). Here we describe the selection procedure of objects for spectroscopic ob- 
servations, the observations themselves, the data reduction, and the construction of the spec- 
troscopic catalogue containing roughly 500 galaxies with secure redshifts. Furthermore we 
discuss global properties of the sample like its distribution in colour-redshift space, the ac- 
curacy of redshift determination, and the completeness function of the data. We derive the 
rest-frame if' -band luminosity function of galaxies at median redshifts of z = 0.2, z = 0.4, 
and z = 0.7. We find evidence for mild evolution of magnitudes (AM^- ~ —0.70 mag) and 
number densities ( / <fr* K ~ —0.35) to redshift one. Furthermore, we present the rest- 
frame J-band luminosity function of galaxies at these redshifts, the first determination of this 
quantity at higher redshifts, with a behaviour similar to the if -band luminosity function. 

Key words: surveys - infrared: galaxies - galaxies: photometry - galaxies: evolution - galax- 
ies: luminosity function - galaxies: fundamental parameters - cosmology: observations 



1 INTRODUCTION 

The study of the luminosity function of galaxies at different wave- 
lengths and at different cosmic epochs is one of the most important 
methods to address the problem of the formation and evolution of 
galaxies within observational cosmology. 

During the last decades, most field galaxy surve ys have be en 
optically selected, mostly in the I or B band (see e.g. Ellis 1997 for 
a review). In recent years however, there has been growing interest 
in the study of near-infrared-selected samples of galaxies. Espe- 
cially the it -band at wavelengths of roughly A ~ offers the 
unique opportunity to detect evolved galaxies, since the if -band 
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light of galaxies is largely dominated by radiation from the evolved 
stellar population iRix & Riekell993tlKauffmann & Charloll99l 
Brinchmann & Ellis 2000). Moreover, in the case of the if -band, 
the transformation which has to be applied to convert the observed 
flux of an object's redshifted spectral energy distribution through a 
given filter to the corresponding rest-frame flux, the k(z) correc- 
tion, is much smaller - and indeed negative - than in the optical 
wavebands, and the variation of k(z) corrections between differ- 
ent galaxy types is small. These effects reduce systematic errors 
introduced by applying k(z) corrections during the construction 
of the luminosity function considerably. As another advantage of 
near-infrared selected surveys, the rest-frame A'-band luminosity 
can be used to deriv e stellar masses for the galaxies, as has been 
done bv lDrorv et all (2001b, hereafter MUNICS III), for example. 
Obviously, this allows the direct study of the history of the assem- 
bly of stellar mass in the galaxies and thus helps to understand the 
mechanism of their formation. 

The past ten years have seen a number of measurements of 
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near-infrared luminosity functions, all of which are summarised in 
Tables l"i"2l and IT3l The easiest way to determine the A-band lu- 
minosity function of galaxies is by follow-up near-infrared pho- 
tometry of existing, typic ally optically select ed, g alaxy re dshift 
catalogues. Mobasher. Sharp ies & Ellisl <1993h and Szokolv et al. 
( 1998) have chosen this approach, and later on Lovedav ( 200JJ) has 
presented the luminosity function at a median redshift z — 0.051 
from A-band imaging of fcj-sele cted galax ies from th e Stromlo- 
APM Redshift Survey iMaddox et alJll99dlLovedav et alJll99fj) . 
Note that the use of optically-selected samples can lead to a bias 
against red (for example old or dusty) galaxies. 

For near-infrared selected galaxy surveys, most measure- 
ments of the A-band luminosit y function so fa r were based on 
local (z < 0.15) samples. Gar dner et alj fl997) have calculated 
the A'-band luminosity function at a median redshift z — 0.14. 
iKochanek et~al]<200ll) have determined type-dependent luminosity 
functions at median redshift z = 0.023 from 2MASS l Jarrett et al. 
20001) and the CfA2 jGeller & HuchrJl989l) and UZC <Falco etaf 
catalogues. A similar strategy was adopted bv Tcole^Hd* 
(2001) who combined the photometric 2MA SS d ata with the 
2dF G alaxy Redshift Survey iFolkes et alJl999l) . and lBalogh eTail 
1200 ll) who cross-corre late the 2 MASS data with the Las Cam- 
panas Redshift Survey l Shect manet ai]|l996l) . As will be shown 
below, there is good agreement between these determinations of 
the rest-fr ame A'-band luminosity function in the local universe. 
However, Huan g"et"al] 120021) have presented a measurement of 
the local A-band luminosity function from the Hawaii-AAO A- 
band Galaxy Redshift Survey, finding a slightly brighter Al* and 
a steeper faint-end slope, which they attribute to the influence of 
different redshift ranges of the local samples. They also note that 
their value for the faint-end slope is in better agreement with optical 
luminosity functions. 

At higher redshifts, Glaze brook et alj ll995l) derived the A'- 
band luminosity function out to redshifts z < 0.8 and find ev- 
idenc e for a brig htening of the characteristic luminosity at z > 
0.5. ICowie et alj 1 19961) present the evolution of the A-band lu- 
minosity function of galaxies in four redshift bins z G [0,0.1], 
2 G [0.1, 0.2], z G [0.2, 0.6], and z G [0.6, 1.0] based on a deep, 
but rather small spectrosco pic sample, co mplemented by the shal- 
lower samples from Songai la et alj J 1994) . They find no evolution 
with redshift. Using our own, much larger sample of spectroscop- 
ically calibrated photometric redshifts from the MUNICS survey, 
we find mild evolution to redshift one, with a brightening of 0.5 
to 0.7 mag, and a decrease in number density of about 25 per cent 
iDrorv et alj2"003l hereafter MUNICS II). 

The spectroscopic sample of A-band selected galaxies de- 
scribed in this paper enables us to derive the rest-frame A-band 
luminosity function of galaxies at redshifts z = 0.2, z — 0.4, and 
2 = 0.7 from data bas ed on a survey mu ch larger in area than 
the Hawaii Deep Fields ICowie et alj |1996). making the luminos- 
ity function presented in this paper much less affected by cosmic 
variance. 

Furthermore, the rest-frame J-band luminosity function at 
these redshifts is presented for the first time. So far, there are only 
two local measurements of the luminosity function in this band 
iBaloeh et all200ll : ICole et all200ll) . 

The paper is organised as follows. Section|2|briefly describes 
the Munich Near-Infrared Cluster Survey (MUNICS), the selec- 
tion of objects for spectroscopy, the observations, as well as the 
data reduction and the construction of the redshift catalogue. The 
properties of the spectroscopic sample are discussed in Section[3] 
Section |4| describes the A-band and J-band luminosity functions 



Table 1. The five MUNICS survey fields for which spectroscopic data 
are available. The table gives the field name, the field coordinates for the 
equinox 2000, the limiting magnitudes for the spectroscopic observations, 
and the number of objects observed spectroscopically in each field. The 
field named S2F1 contains the sparse sample observed with the ESO VLT, 
the data in the field S2F5 are rather incomplete and are excluded from any 
further analysis. 



Field a (2000.0) 5(2000.0) Lim. mag. Spectra 



S2F1 


03:06:41 


+00:01:12 


K < 19.0 


347 


S2F5 


03:06:41 


-00:13:30 


K < 17.5 


29 


S5F1 


10:24:01 


+39:46:37 


K < 17.5 


121 


S6F5 


11:57:56 


+65:35:55 


K < 17.5 


193 


S7F5 


13:34:44 


+ 16:51:44 


K < 17.5 


140 



of galaxies as derived from the spectroscopic data, with a detailed 
discussion of the results in Section|5| Finally, Section|6|concludes 
this work with a summary of our results. 

We assume Qm = 0.3, Qa = 0.7 throughout this paper. We 
write Hubble's Constant as Ho — 100 h kms -1 Mpc -1 , using 
h — 0.60 unless explicit dependence on h is given. For conve- 
nience, we write the apparent A'-band magnitude as ttik in some 
places. 



2 SPECTROSCOPIC OBSERVATIONS 

2.1 The Munich Near-Infrared Cluster Survey 

The Munich Near-Infrared Cluster Survey (MUNICS) is a wide- 
field medium-deep survey in the near-infrared A' and J pass- 
bands which is fully described in D rorv et all (2001a; hereafter 
MUNICS I). In brief, the survey consists of a A'-selected cata- 
logue down to K' < 19.5 (50 per cent completeness limit for 
point sources) covering an area of 1 square degree. Additionally, 
0.35 square degrees have been observed in I, R, and V. B-band 
imaging has been recently completed for the same area and will be 
described elsewhere. The limiting magnitudes of this main part of 
the survey are 23.5 in V, 23.5 in R, 22.5 in 7, and 21.5 in J. The 
magnitudes are in the Vega system and refer to 50% completeness 
for point sources. The fields used in this work are five randomly 
selected patches of sky (see TableQfor details). The layout of the 
survey, the observations and data reduction are described in MU- 
NICS I. 

2.2 Selection of the spectroscopic sample 

Objects for spectroscopic observations were chosen from the K- 
band selected photometric catalogue of MUNICS in five survey 
fields, the details of which can be found in Table[j](see MUNICS I 
for a description of the field nomenclature). 

Object selection for spectroscopy was based on two criteria. 
Firstly, we aimed at a A-band magnitude-limited sample. Due to 
the use of optical spectrographs, the appropriate A-band limit is de- 
termined by the typical colours of red galaxies (roughly R—K ~ 4, 
see Fig.|5J and the limits of the optical spectrographs at the tele- 
scopes we used. Trying to keep the A-band completeness of the 
spectroscopic sample as high as possible yields sample limits of 
A < 17.5 for spectroscopic observations at the Calar Alto 3.5-m 
telescope and A < 19.0 for observations at the Very Large Tele- 
scope (VLT). Obviously, a small fraction of very red objects will be 
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Table 2. Summary of results of the test of the image-based classification 
method for a sample of objects with K' < 16.5 in two survey fields (upper 
part of table) and for the whole spectroscopic sample (lower part of table). 
The numbers quoted in the table give the number and fraction of objects 
with stellar or galactic spectrum among objects classified as point-like or 
extended, respectively. 



Morphological classification 




Spectral classification 






Star 


Galaxy 


Sparse sample with K' < 16.5 








Point-like 


39 


100% 


0% 


Extended 


11 


20% 


44 80% 


Complete spectroscopic sample 








Point-like 


64 


79% 


17 21% 


Extended 


28 


7% 


397 93% 



lost, but thei r number densi ty is comparatively small anyway (see, 
for instance, Martini (2001) and references therein). Results on the 
few very red objects identified spectroscopically can be found in 
Sectionl3~8l 

Secondly, in selecting objects for spectroscopy, we have tried 
to exclude stars. This was done using the image-based classification 
of objects into point-like objects and extended sources as described 
in MUNICS I. 

To test this classification procedure, we have built up a small 
spectroscopic test sample of bright objects which was purely mag- 
nitude selected and therefore contains both stars and galaxies. The 
results of the comparison between image-based and spectroscopic 
classification show that our morphological approach is able to dis- 
tinguish stars and galaxies with reasonable reliability. The results 
of this test are summarised in Table|2| 

On a purely magnitude-limited sample containing all objects 
with K' < 16.5, all objects classified as point-like are indeed stars, 
and the stellar contamination is 20 per cent only. On the other hand, 
we can also use the complete spectroscopic catalogue to investigate 
the reliability of our classification method, yielding also very good 
agreement between spectral and morphological classification. 

Thus our method to classify objects into point-like or extended 
sources can reliably reduce stellar contamination of the spectro- 
scopic catalogue as long as the objects are not too faint. Pre- 
selecting our spectroscopic sample by image morphology results 
in a loss of less than 4 per cent of all galaxies which are classified 
as point-like, and in a 7 per cent contamination by stars. 

2.3 Spectroscopic observations and data reduction 

The largest part of the spectroscopic observations was carried out 
with the Multi-Object Spectrograph for Calar Alto (MOSCA) at 
the 3.5-m telescope at Calar Alto Observatory (Spain), and with 
FOcal Reduc er and low-dispersion Spectrograph (FORS) 1 and 2 
<Seifertetal]ll994l) at the European Southern Observatory's Very 
Large Telescope (VLT). Part of the sample was obs erved with the 
Low Resolution Spectrograph (LRS: lHill et alll998i) at the Hobby- 
Eberly Telescope (HET) at McDonald Observatory, Texas, and with 
ESO Faint Object Spectrograph and Camera (EFOSC) 2 at the ESO 
3.6-m telescope on La Silla (Chile). All observing runs are listed in 
Table[3] 

MOSCA was used with the Green 500 grism and without 
any filter. MOSCA's multi-object spectroscopy mode uses slit 



Table 3. Observing runs of spectroscopic follow-up observations for the 
MUNICS project. The runs in 2002 were mostly carried out in service 
mode. 



Date 


Telescope 


Instrument 


15.12.1999 


ESO 3.6 


EFOSC2 


26.-31.5.2000 


Calar Alto 3.5 


MOSCA 


27.5.2000 


Hobby-Eberly Telescope 


LRS 


21.11.2000 


ESOVLTUT1 (Antu) 


FORS1 


21.-22.11.2000 


ESO VLT UT2 (Kueyen) 


FORS2 


24.-28.11.2000 


Calar Alto 3.5 


MOSCA 


17.-20.1.2001 


Calar Alto 3.5 


MOSCA 


26.3.-1.4..2001 


Calar Alto 3.5 


MOSCA 


18.-21.5.2001 


Calar Alto 3.5 


MOSCA 


15.-20.12.2001 


Calar Alto 3.5 


MOSCA 


9.-12.4.2002 


Calar Alto 3.5 


MOSCA 


9.-17.5.2002 


Calar Alto 3.5 


MOSCA 


7.10.2002 


Calar Alto 3.5 


MOSCA 


9.-10.11.2002 


Calar Alto 3.5 


MOSCA 



Table 4. Technical parameters of the spectrographs used for the observa- 
tions. The effective resolution of the resulting spectra is given for the grisms 
listed in the table and the appropriate slit-widths mentioned in the text (the 
seeing during the spectroscopic observations always was of the order of the 
slit- width). 



Instrument 


Grism d 


k Filter 


Spectral range 


Resolution 


MOSCA 


Green 500 


4300 A - 8000 A 


13.6A 


EFOSC2 


Gil 




3380 A -7520 A 


13.2 A 


LRS 


G300 4 


- GG385 


4000 A - 8000 A 


13.9 A 


FORS1 


G300I - 


f OG590 


6000 A - 9500 A 


13.0A 


FORS2 


G300I - 


f OG590 


6000 A - 9500 A 


13.0A 



masks containing typically 20. . .25 slits of 1.5" width. MOSCA 
is equipped with a 2048x4096 pixel CCD with 15/^m pixel size, 
yielding an effective area of 10' x 10' usable for spectroscopy. 

At the VLT, FORS1 offers 19 movable slit-lets, whereas 
FORS2 is equipped with a Mask Exchange Unit (MXU), which 
allows spectroscopy of a larger number of objects. For both in- 
struments a slit width of 1" was chosen. The CCD detectors are 
2048 x 2048 pixel in size with 24/xm pixels, thus allowing a field of 
roughly 3'x 7' to be used for spectroscopy. Grism 300 I and filter 
OG590 was used for the observations 

In multi-object spectroscopy mode, the LRS at HET provides 
13 slit-lets of 1.5" width. The spectra were obtained through grism 
300 and filter GG385. 

Finally, EFOSC2 at the ESO 3.6-m telescope was equipped 
with grism 11 and a slit- width of 1.0 ". The technical characteris- 
tics of all instruments used for the spectroscopic observations are 
summarised in Table|4] 

The objects selected for spectroscopy were sorted into bins 
according to their apparent i?-band magnitude. The magnitude 
ranges, the typical number of mask setups for each MUNICS field 
(13' x 13'), and the exposure times are given in Table [5] All ob- 
jects with R < 21.5 (roughly corresponding to K' < 17.5) were 
observed at the Calar Alto 3.5-m telescope, the ESO 3.6-m tele- 
scope, and the HET, whereas spectroscopy of the fainter sample 
was carried out at the ESO VLT. 
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Table 5. i?,-band magnitude bins, typical number of masks per MUNICS 
field, and exposure times for the spectroscopic observations. Note that the 
faintest objects with R > 21.5 were observed at the ESO VLT, while all 
brighter objects were mainly observed at the 3.5-m telescope at Calar Alto, 
a few also at the ESO 3.6-m telescope, or at the HET. 



Magnitude 


Number of masks 


Exposure time 


R < 18.5 


1 


1 X 1800 s 


18.5 < R < 19.5 


2 


2 X 2100 s 


19.5 < R < 20.5 


3 


5 X 2400 s 


19.5 < R < 21.5 


4 


5 x 5400 s 


R > 21.5 


9 


3 x 3000 s 



Data reduction was performed using IRAF , except for cosmic 
ray filtering. Firstly, a standard CCD reduction was performed, in- 
cluding over-scan and bias correction, as well as interpolation of 
bad pixels. Cosmic rays were identified by searching for narrow 
local maxima in the image and fitting a bivariate rotated Gaussian 
to each maximum. A locally deviant pixel is then replaced by the 
mean value of the surrounding pixels if the Gaussian obeys appro- 
priate flux ratio and sharpness criteria <G6ssl & Riffesej2002l) . 

The following reduction of the spectra was performed within 
the IRAF SPECRED package. The locations of the spectra on the 
two-dimensional image were defined, and the trace of the spectra 
produced by the distortion of the spectrograph was approximated 
by a polynomial of fourth order. Flat-fields were constructed from 
typically 21 images of an internal continuum lamp, and normalised 
by fitting a set of low-order splines to the illumination function 
of each slit. After dividing the frames by the appropriate flat-field, 
the apertures for extracting one-dimensional spectra from the two- 
dimensional images as well as the apertures for subtraction of the 
night sky were defined. This was done manually for each object 
to ensure maximum signal-to-noise ratio and accurate sky subtrac- 
tion. If possible, sky apertures on both sides of each object were 
used, and the variation of the night sky perpendicular to the spec- 
trum of each object was approximated by a polynomial of second 
order. Note that we used lower-order polynomials in cases where 
there were only few pixels to define the sky region in order to pre- 
vent artefacts due to ill-defined approximations. After that the one- 
dimensional spectra were extracted by summing the pixel values for 
each wavelength bin within the manually defined apertures and sub- 
tracting the appropriately scaled sky value from the fit perpendic- 
ular to the spectrum. Additionally we produced two-dimensional 
sky-subtracted images to be able to check whether a feature seen in 
the one-dimensional spectrum is a real spectral feature or a resid- 
ual, for example from the replacement of cosmic ray events or the 
subtraction of sky-lines. 

In the case of the MOSCA, EFOSC2 and LRS observations, 
spectra of internal lamps were used for wavelength calibration. The 
wavelength transformation was computed in form of a fourth order 
polynomial. Then the one-dimensional spectra were transformed to 
wavelength coordinates and re-binned using the appropriate reso- 
lution. Note that we usually took only one wavelength calibration 
image per mask, either between two mask exposures or after finish- 
ing all mask exposures. Nevertheless, this does not severely affect 

1 IRAF, the IMAGE REDUCTION AND ANALYSIS FACILITY, is distributed 
by the National Optical Astronomy Observatories, which are operated by 
the Association of Universities for Research in Astronomy, Inc., under co- 
operative agreement with the National Science Foundation. 



the wavelength calibration, since the maximum shift of the spectral 
pattern of a calibration lamp, caused by the flexure of the instru- 
ment, was measured to be 1 pixel for MOSCA between two ex- 
treme orientations of the spectrograph. According to the manuals, 
the flexure of the other instruments should have a similar order of 
magnitude. Also, the accuracy of the wavelength calibration can be 
efficiently examined by measuring the positions of night-sky emis- 
sion lines. The comparison showed negligible deviations. The influ- 
ence on the redshift determination will be discussed in Section l7l4l 
For the FORS data we used the lines of the night sky for the wave- 
length calibr ation, with the line data taken from Osterbrock et al. 
h99flll997l) . 

Finally, the spectra were flux calibrated using spectro- 
photometric standard stars observed during each night. Extinction 
correction was based on average extinction curves either as pub- 
lished by the observatories (ESO for La Silla and Paranal Observa- 
tory, and McDonald Observatory) or by Hope & Fernandez 1 2002) 
for the Calar Alto observations. The accuracy of the flux calibration 
was tested by comparing the flux in individual broad-band photom- 
etry filters to the flux in the flux-calibrated spectrum. Fig.Qshows 
an example for this method. Clearly, the flux calibration is satisfy- 
ing, with deviations visible only at the ends of the spectral range of 
the instrument, where the corrections due to the response function 
of the grism become large. 

Redshifts were determined using the RV package within IRAF. 
Prominent absorption and emission features were identified in the 
spectra to obtain the corresponding redshift of a galaxy. 2 We com- 
pared the one-dimensional spectrum with the two-dimensional sky- 
subtracted image to exclude features possibly affected by residuals 
from night-sky subtraction or cosmic filtering. Furthermore we en- 
sured that the radial velocity displacement between different spec- 
tral features is not larger than the typical error expected from the 
measurement of the line centre and the wavelength calibration, i.e. 
typically Az <5 x 10" 4 . 



2.4 Catalogue construction 

The final spectroscopic catalogue contains the celestial coordinates, 
the pixel coordinates form the photometric catalogue, the instru- 
mental identification (mask and slit number), the measured redshift, 
the object classification based on the observed spectrum (galaxy, 
AGN, star, or unidentified object), a confidence class for the red- 
shift determination, and a list of spectral features used for mea- 
suring the object's redshift. The confidence class assigned to each 
object ranges from "1" to "6", where "1" means highest confi- 
dence, and "6" means that no redshift could be determined. De- 
scriptions of the confidence classes can be found in Table [6] al- 
though these can give only a crude impression of the meaning of the 
confidence classes. For all further analysis objects with confidence 
classes "5" and "6" are excluded. To check the assignment of confi- 
dence classes, we have made use of the 44 repeat observations of 41 
objects described in detail in Section l7l4l The agreement between 
confidence classes attributed to the same object from different ob- 
servations is in general very good. Deviations can be explained by 



2 So far we have not applied cross-correlation techniques to determine red- 
shifts, however we plan to do so in the future to reduce the fraction of galax- 
ies without secure redshift due to low signal-to-noise ratios, although first 
tests show that we cannot significantly increase the number of securely mea- 
sured redshifts by using correlation methods. 
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Figure 1. Comparison of the flux calibrated spectrum (solid line), the flux derived from MUNICS broad-band photometry in the filters B, V, R, I, J, and 
K (filled circles with error bars), and a fit of a stellar population synthesis model to the photometric data (dotted line) for two objects in the spectroscopic 
MUNICS catalogue. 



Table 6. Description of the confidence classes for the redshift measurement 
used in the final spectroscopic catalogue. 



Class Description 

1 Many spectral features securely identified, very high 
signal-to-noise ratio 

2 Many spectral features securely identified, high 
signal-to-noise ratio 

3 Many spectral features securely identified, intermedi- 
ate signal-to-noise ratio 

4 Several spectral features identified, lower signal-to- 
noise ratio 

5 Tentative identification of only a few lines or objects 
with single emission line, low signal-to-noise ratio 

6 no redshift determination possible (either very low 
signal-to-noise ratio or technical problem) 



3 PROPERTIES OF THE SPECTROSCOPIC SAMPLE 

3.1 Object classes 

The final catalogue contains 830 objects. 725 (or 87.3 per cent) are 
galaxies, 3 (0.4 per cent) are active galactic nuclei, 100 (12.0 per 
cent) are stars, and 2 objects (0.2 per cent) remain unidentified due 
to bad quality of the spectra or a lack of prominent spectral features 
in the wavelength range covered by the spectra. Redshifts could be 
determined for 500 objects (confidence classes 1 to 4), or 68.7 per 
cent of all extragalactic objects. 

3.2 Confidence classes 

The number of objects in the confidence classes described above 
are given in TableQ The first three confidence classes contain 57.1 
per cent of all objects selected for spectroscopy; these objects have 
high-quality spectra. Redshifts could be determined for 68.7 per 
cent of the objects. 



different signal-to-noise ratios due to weather conditions or tech- 
nical reasons. Not a single object with confidence class "1" to "4" 
has got a deviant redshift from a repeat observation. Typical exam- 
ples for spectra and the appropriate confidence classes are shown 
in Fig.|2| 



3.3 Redshift distribution 

The redshift distribution of the spectroscopic sample is shown in 
Fig. [3] The histogram has the shape expected for a magnitude lim- 
ited spectroscopic sample, except for the slight excess of objects 
with redshifts z > 0.6 which, of course, are due to the sparse sam- 
ple of faint objects observed at the VLT. The lack of objects with 
redshifts from VLT spectroscopy in the range 0.5 < z < 0.6 is 
most likely due to the limits of the instruments' spectral range. Ob- 
jects with these redshifts have their 4000A break around the lower 
end of the spectral range and the Ha line around the upper end 
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Confidence class 1 [2301_01] 



Confidence class 2 [6403 01] 



Confidence class 3 [5401_04] 




4500 5000 5500 6000 6500 7000 7500 8000 
X [Angstrom] 




4500 5000 5500 6000 6500 7000 7500 8000 
X [Angstrom] 




4500 5000 5500 6000 6500 7000 7500 8000 
X [Angstrom] 



Confidence class 4 [210109] 



Confidence class 5 [6302 1 3] 



Confidence class 6 [5401_23] 
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Figure 2. Examples for the assignment of confidence classes, ranging from "1" (upper left) to "6" (lower right). The numbers in square brackets denote the 
internal identification number of the object. 



Table 7. Distribution of extragalactic objects into the 6 confidence classes 
used in the final spectroscopic catalogue. 



Class 


Number of objects 


Per cent 


1 


162 


22.25 


2 


137 


18.82 


3 


117 


16.07 


4 


84 


11.54 


5 


8 


1.10 


6 


220 


30.22 


Total 


728 


100.00 



of the spectral range. Thus redshifts of objects with low signal-to- 
noise ratio may be difficult to determine due to the absence of those 
prominent spectral features in the observed spectra. This kind of 
bias can be corrected for by looking at the fraction of galaxies with 
securely determined redshift as a function of two colours, as will 
be shown below. 
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Figure 3. The overall redshift distribution of all galaxies and active galactic 
nuclei in the final spectroscopic catalogue (solid line), divided into objects 
with spectroscopy from the Calar Alto 3.5-m telescope, the ESO 3.6-m tele- 
scope or the HET (dashed line) and from the VLT (dotted line). Note the 
excess of objects at redshifts above 0.6 due to the VLT observations. 



3.4 Accuracy of redshift determination 

For objects with secure redshifts, where we consider all spectral 
features to be correctly identified, possible sources of error in the 
determination of redshift are firstly a possible displacement be- 
tween the object and the slit centre, secondly the error of the wave- 
length calibration, and finally the error in the measurement of the 
line centres. 

(i) An offset of the object with respect to the slit centre causes an 
error of the redshift measurement. We estimate the RMS error to be 
0.5", corresponding to an inaccuracy of a ~ 6 x 1CP 4 in redshift 
space. This value is estimated from the accuracy of astrometry in 
the MUNICS catalogue (see MUNICS I), the typical errors of mask 
alignment and under the reasonable assumption that errors in the 



mechanical production of the masks (or in the slit-let positioning 
for FORS1 and the LRS) can be neglected. 

(ii) The random errors of the wavelength calibration, the RMS 
of the residuals of the fit to the identified lines, are of the order 
0.3 A, or, in redshift, a ~ 4 x 10~ 5 . Compared to the other sources 
of error, this inaccuracy can be neglected. For MOSCA, the maxi- 
mum shift of the spectrum of the calibration lamp between zenith 
and airmass 2 was measured to be 1 pixel, corresponding to a sys- 
tematic error of a < 4 x 10~ 4 in redshift. 

(iii) Typical errors of the measurement of the line centres for the 
redshift determination are <r < 5 x 10~ 4 . 

The overall limit to the accuracy estimated from these errors 
is a ~ 8 x 10~ 4 in redshift space. 
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Figure 4. Histogram of the absolute values of pairwise differences of the 
measured redshifts of objects with multiple spectroscopic observations in 
the catalogue (solid line). The dotted line shows a Gaussian approximation 
to the histogram for which a' = 12 X 10 — 4 , corresponding to a an RMS 
eiTor of single measurements of<r = 8xl0 -4 . 

We can also empirically determine the error of the redshift 
measurements from multiple observation of objects. A total of 41 
galaxies - about 8 per cent of the sample - has been observed re- 
peatedly, usually in order to fill gaps on the slit masks used for 
observation, thus enabling us to estimate typical errors of the red- 
shift measurement by comparing the redshifts obtained from differ- 
ent observations. Note that 13 of these objects have been observed 
both with MOSCA and FORS, and 5 using MOSCA and the LRS. 
Thus it is possible to check for any systematics between redshift 
determinations with different spectrographs, and, in particular, for 
systematic offsets between different means of wavelength calibra- 
tion, since we used the lines of the night sky for the FORS obser- 
vations, whereas spectra from calibration lamps were used for the 
other spectrographs. 

The histogram of the absolute values of pairwise differences 
for the 44 repeated redshift measurements for the 41 objects is 
shown in Fig. |4| Also shown is a Gaussian approximation to the 
histogram having a width characterised by a' = 12 x 1CU 4 , cor- 
responding to a an RMS error of single measurements of a = 
<t'/V2 = 8 x 1(T 4 , exactly the same value we have derived 
from the formal error analysis. Furthermore, we do not find any 
systematical differences between redshift determinations achieved 
with data from different instruments or observing runs. 

3.5 Colours, magnitudes, and redshifts 

To give a feeling for the distribution of objects in magnitudes, 
colours, and redshifts, we present diagrams for various combina- 
tions of these properties in Fig. [5] 

In the upper panels of the figure, we show distributions of the 
apparent if-band magnitudes and the R — K colours versus red- 
shift. The comparison of model spectral energy distributions from 
a combina tion of empirical spectra and stellar population synthesis 
models by Maraston ( 1998) for various Hubble types to the R—K 
colours of the objects in the MUNICS spectroscopic catalogue (up- 
per right panel of Fig.|5J shows good agreement, although we seem 
to miss blue objects at higher redshifts, a fact which must be ac- 
counted for in any analysis of this dataset. 

The two lower panels of Fig. [5] shows the distribution of 
extended objects in the photometric catalogue together with the 



colour distribution of galaxies with spectroscopic redshifts. Com- 
parison shows that the spectroscopic catalogue is a fair represen- 
tation of the distribution of objects in the photometric sample. 
Clearly, there is some incompleteness for faint red objects, but this 
can be corrected for, as will be shown in the following Section. 



3.6 Redshift sampling rate and sky coverage 

Observations for the spectroscopy of objects from the photomet- 
ric MUNICS catalogue are almost complete, and most of the fields 
have been observed with reasonable completeness. The fraction of 
galaxies with redshift among all galaxies in the photo metric cat- 
alogu e is usually called the redshift sampling rate (e.g. iLin et all 
1999). As an example, Fig.|6|shows the redshift sampling rate of 
objects in all survey fields with spectroscopic data. The redshift 
sampling rate is the fraction of objects with successful redshift de- 
termination among all galaxies in the photometric catalogue, and, 
of course, is different from the redshift success rate, which is the 
fraction of all galaxies with redshift among all galaxies in the spec- 
troscopic catalogue. The redshift success rate depending on appar- 
ent magnitude and colours is shown in Fig.Qand discussed in Sec- 
tion |3J| These colour distribution should be compared to the dis- 
tributions of objects from the photometric and spectroscopic cata- 
logue shown in Fig. [5] 

Although a few of the masks for the fields still have to be ob- 
served, this does not necessarily inhibit the use of the sample for 
analysis of, e.g., the luminosity function, as long as the sample of 
objects with successful redshift determination is a fair represen- 
tation of the total sample, and as long as any systematic incom- 
pleteness effects can be corrected for. We will show that this is the 
case for our catalogue because of the procedure we followed during 
mask preparation and observation. Firstly, we have tried to select 
all objects with K' < 17.5 in five survey fields for spectroscopy at 
Calar Alto. Since the field of view of the spectrograph (MOSCA) 
is almost as large as the size of our survey field, and since we usu- 
ally have several masks per magnitude bin (at least for the fainter 
objects), we are not limited by geometrical constraints from the ar- 
rangement of slits on each mask. For the brighter magnitude bins 
we sometimes had to drop a few objects during the preparation of 
a mask, but we have tried to include these objects in fainter mask 
setups. Hence we expect that the small fraction of objects on to 
which no slit could be put is statistically similar to the distribution 
of objects to the appropriate magnitude limit. Secondly, all masks 
which could be observed so far were selected randomly from the 
masks prepared for observations, so that we do not expect any se- 
lection bias here either. Fig.|8|shows the distribution of objects with 
successful spectroscopy as compared to objects in the photomet- 
ric catalogue for the survey patches called S2F1, S5F1, S6F5, and 
S7F5. These particular fields have highest completeness. The field 
S2F1 contains the sparse sample of fainter objects observed with 
the VLT. The figure nicely illustrates the statistically homogeneous 
distribution of the objects in the spectroscopic catalogue. The ex- 
cess of objects with successful redshift determination visible in the 
lower half of the field S2F1 is due to foreground structure which 
can readily be seen in the image itself, showing again the influence 
of cosmic variance and the importance of having several fields with 
a relatively large solid angle for the analysis of statistical quanti- 
ties, one of the huge advantages of the MUNICS project compared 
to smaller surveys. 
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Figure 5. Upper left panel: The distribution of apparent ii"-band magnitudes of galaxies versus redshift in the MUNICS spectroscopic sample. Upper right 
panel: The R— /("-colour versus redshift diagram for galaxies with spectroscopic redshifts. Also shown are model spectral energy distributions from (see text 
for details), roughly corresponding to Hubble types E (solid line), Sa (dash-dotted line), Sb (dashed line), and Sc (dotted line). Lower left panel: The J — K 
versus K colour-magnitude diagram for galaxies (filled circles) and stars (filled triangles) with good-quality spectra. The open circles show the distribution 
of objects classified as extended in the photometric catalogue of one of the MUNICS fields (S2F1). Also indicated is the magnitude limit of the large and 
bright spectroscopic sample at K = 17.5 (dotted line). Lower right panel: Same as lower left panel, but for the distribution of R—K colours versus apparent 
A"-band magnitude. 
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Figure 6. The redshift sampling rate, i.e. the fraction of successful redshift determinations among galaxies in the photometric catalogue, as a function of m K 
(left panel), R—K (middle panel), and J—K (right panel) for all galaxies in the survey patches S2F1, S5F1, S6F5, and S7F5 (the field S2F5 has been excluded 
from further analysis because of the small number of available spectra). The dotted line in the left panel indicates the formal limit of K < 17.5 of the main 
part of the spectroscopic survey. The colour distributions in the middle and the right-hand panel are those of objects brighter than this limit. 
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Figure 7. The redshift success rate, i.e. fraction of successful redshift determinations among all galaxies which were observed spectroscopically, as a function 
of mji (without any colour constraint; left panel), R—K (middle panel), and J—K (right panel) for all galaxies in the survey patches S2F1, S5F1, S6F5, and 
S7F5). The dotted line in the left panel indicates the formal limit of K < 17.5 of the main part of the spectroscopic survey. The colour distributions in the 
middle and the right-hand panel are those of objects brighter than this limit. 



3.7 Redshift success rate 

The efficiency of redshift determination of spectroscopic observa- 
tions is described by the redshift success rate rather than the red- 
shift sampling rate. The former is the fraction of objects with secure 
redshift among all galaxies which were observed spectroscopically 
and is shown in Fig.0 As expected, the redshift success rate drops 
off for very faint and very red objects, but is in general very high 
due to the good quality of the spectra. 

3.8 Extremely red objects 

During the last years there has been a lot of research on 'extremely 
red objects (EROs)', usually defined in terms of their very red 
optical -near-infrared colour {R — K > (5 ... 6) mag; see, for in- 
stance, iMartinilEoOll) and references therein). While the MUNICS 
catalogue obviously contains such objects, they are not aimed at 
with the spectroscopic observations described in this paper, mainly 
because it is very difficult to obtain optical spectra of faint EROs. 
This is evident from Figures[6]and[7]showing the redshift sampling 
rate and the redshift success rate for the MUNICS spectroscopic 
observations. Nevertheless, there are 15 such objects in the spec- 
troscopic catalogue, all having colours of 5.0 mag < R— K < 
5.5 mag. Among these we find 7 objects for which no redshift 
could be determined, 5 galaxies with redshifts 0.46 < z < 1.01, 
all having a spectral energy distribution characteristic of early-type 
galaxies, as well as 3 stars of spectral type M. First results of near- 
infrared spectroscopy of EROs selected from the MUNICS cata- 
logue are described in lSaracco et alj f2003). 



4 THE NEAR-INFRARED LUMINOSITY FUNCTION OF 
GALAXIES 

4.1 Calculation of the luminosity function 

The luminosity funct ion is com puted using the non-parametric 
Knax formalism l Schmidt 1968). This method has been shown to 
yield an unbiased estimate of the luminosit y function if t he sample 
is not affected by strong clustering l Takeuc hi et all200oT) . Because 
of our field selection and the relatively large area of the survey, di- 
vided into several individual fields, we believe that this assumption 
is valid for our sample. 



The Vmax formalism accounts for the fact that some fainter 
galaxies are not visible in the whole survey volume. Each galaxy 
i in a given redshift bin [ziowor , z U p P cr] contributes to the number 
density an amount inversely proportional to the volume Vi in which 
the galaxy is detectable in the survey: 

min(z uppcr ,z ma x) 

Vi = J ^dz, (1) 

z lowcr 

where dV = dil r 2 dr is the comoving volume element, dQ, 
is the solid angle covered by the survey, and z m ax is the maximum 
redshift at which galaxy i having absolute magnitude Mi is still 
detectable given the limiting apparent magnitude of the survey. We 
have made sure that the effect of the volume correction is of impor- 
tance only in the faintest bin in absolute magnitude, and that even 
in this case the correction is small. 

Additionally, the contribution of each galaxy i is weighted by 
the inverse of the detection probability P(mic,i) of the .ftf'-band se- 
lected photometric catalogue, where we assume that the detection 
probability is independent of the galaxy type and can be approxi- 
mated by that of point-like sources. Completeness simulations for 
realistic galaxy profiles at various redshifts have shown that this 
approximatio n is ind eed sufficient for galaxies at redshifts z < 1 
I Snigul a et ai]|200l hereafter MUNICS IV). However, since the 
objects under consideration here are comparatively bright, the in- 
fluence of this correction is negligible anyway. 

In addition to the correction for the incompleteness of the pho- 
tometric MUNICS catalogue described above, we have to correct 
for the incompleteness of the spectroscopic catalogue with respect 
to the photometric sample, described by the redshift sampling rate 
(see Section l3~6t . In principle, this correction depends on the ap- 
parent magnitude of the objects (faint objects may produce a lower 
signal-to-noise ratio or might not have been considered for spec- 
troscopy) 3 , on their intrinsic type (it is easier to determine a redshift 
for objects showing prominent emission lines, for example) and on 
the redshift of the source (influencing the position of prominent 

3 Note that we can neglect any influence of the size of the objects (large 
objects might have larger losses of light at the slit of the spectrograph), 
since - especially at faint magnitudes - objects appear to be almost point- 
like. 
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Figure 8. Sky coverage of spectroscopy in the MUNICS fields S2F1 (upper left panel), S5F1 (upper right panel), S6F5 (lower left panel), and S7F5 (lower right 
panel). The small open circles are objects in the photometric catalogue, the distribution of which is governed by the MUNICS field geometry (see MUNICS I). 
All galaxies with K' < 17.5 and successful redshift determination are marked by filled circles; they seem to be uniformly distributed on the sky. Note that 
spectroscopy of the galaxies in the field S7F5 is not complete yet. Coordinates are pixel coordinates in the A"'-band MUNICS frame. 



spectral features with respect to the spectral range or bright night- 
sky emission lines, for example). However, it is difficult to deter- 
mine a completeness ratio depending on spectral type and redshift, 
because this information is lacking for all objects without secure 
redshift measurement. 

Hence the redshift sampling rate is often quantified in terms 
of its dependence on apparent magnitude and two colours instead 
of spectral type and redshift (see Lin et al. 1 1999), for example). In 
this work we compute the redshift sampling rate d, 

(2) 



a = 



N, 



depending on apparent if'-band magnitude tuk, and the 
colours R—K and J—K. Specifically, for each galaxy i with a red- 
shift, we determine the number Ni of galaxies in the photometric 
sample and the number iV z ,i of galaxies with successful redshift de- 
termination in joint bins of apparent magnitude and colours, where 
we count all galaxies j obeying 

\rriK,i — rnK,i ; l ^ 0.70 mag, 
\(R—K)i — (R—K)j\ < 0.70 mag, (3) 
\(J-K)i - (J-K)j\ < 0.35 mag. 

On the one hand, the size of the intervals should be large 



enough to contain a reasonable number of objects to avoid large 
fluctuations due to small-number statistics. On the other hand, the 
intervals should not be too large in order to be able to follow the 
change of the sampling rate. The numbers given in equation J5J 
have been chosen after careful tests. Note that we have used larger 
intervals for uik and R — K since these distributions are broader 
than the one for J—K. To illustrate the general behaviour of the red- 
shift sampling rate, we show projections of the function in Fig. [6] 
Note that we have to apply a correction factor of 0.93 to the sam- 
pling rate because of stellar contamination of the class of extended 
objects, a bias introduced by our image-based object classification 
algorithm, described in detail in Section lT2l 

The near-infrared luminosity function $(M) is then com- 
puted according to the formula 

$(M) dM = V 1 — ±— JL d M, (4) 
^—^ Vi R{m K ,i) Ci 



where the sum runs over all objects i in the redshift range for 
which we want to calculate the luminosity function. 

4.2 Converting to absolute magnitudes 

For the conversion of apparent magnitudes m into absolute magni- 
tudes M we need k(z) corrections for the galaxies, defined by 
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M = m - 5 log f ) - 2.5 log (1 + z) - k(z), (5) 

where <1l(z) is the cosmological luminosity distance of the 
galaxy, and z is the redshift. 

The k(z) corrections for the objects in the spectroscopic cat- 
alogue are obtained by fitting the model spectral energy distribu- 
tions shown in the left panel of Fig. [9| to the broad-band photom- 
etry of the objects. The template spectral energy distributions are 
constructed by combining empirical spectra and stellar population 
synthesis models by Maraston 1 1998). The detailed description of 
these models will be found in MUNICS II. The k(z) corrections for 
the K band and the J band are presented in the middle and right 
panel of Fig. [5] respectively, while Fig. \\\ shows a comparison of 
flux-calibrated spectrum, broad-band photometry, and fitted model 
spectral energy distribution for two objects in the spectroscopic cat- 
alogue. 

We have also tested the influence of using only one 
intermediate-type model spectral energy distributions for deriving 
the k(z) corrections. As is obvious from Fig. [5] this should work 
well for the K band, where the spread between k(z) corrections for 
different models is small, but a bit less so well for the J band with 
its larger variations. Indeed, we do see hardly any difference for the 
if-band luminosity function, and only a small difference in the J 
band. However, it is important to note that in the following we do 
not use k(z) corrections from one model, but those from the fitting 
of template spectral energy distributions to the six-filter broad-band 
photometry of the MUNICS catalogue. 



faintest magnitude bin does not change the result of the fit signifi- 
cantly. 

Note that the agreement between our measurement and the lo- 
cal Schechter functions is also very good at the faint end, although 
with poor statistics and large correction factors due to the incom- 
pleteness. From this we draw the conclusion that our method to 
correct for incompleteness yields reliable estimates. 

At redshifts 0.3 < z < 0.6, we find a mild evolution 
of the luminosity function, with a higher characteristic luminos- 
ity and a lower normalisation, yielding Schechter parameters of 
M£-51og/i = (-24.04±0.26) magand$K = (0-71 ±0.25) x 
10 -2 h 3 Mpc -3 . The error contours for the Schechter parameters 
for the luminosity function in the two lower redshift bins are shown 
in the same figure. These contours were computed from the \ 2 dis- 
tribution. Finally, we show a comparison of the measurements of 
the if -band luminosity function in all three redshift bins. Although 
the statistics becomes rather poor in the highest redshift bin with 
median redshift z = 0.7, it confirms the trend for the evolution of 
the luminosity function. However, it is evident that the total evolu- 
tion out to redshifts around 0.7 is not dramatically large. 

C omparison of our result with the measurements by 
ICowie et alJ J 19961) who derived the if -band luminosity function 
at various redshifts < z < 1 from a spectroscopic sample in 
two small fields shows good agreement with our results. Their val- 
ues for the Schechter parameters M^, &*k< an d a K are similar to 
measurements for local samples. Again, our data favour a slightly 
larger characteristic luminosity, but the trend of a falling normali- 
sation with redshift is found both in the MUNICS and in the Cowie 
et al. measurements of the luminosity function. 



4.3 The if -band luminosity function 

The rest-frame if -band luminosity function of galaxies drawn from 
the spectroscopic sample was constructed in the redshift intervals 
0.1 < z < 0.3 (median redshift z = 0.2), 0.3 < z < 0.6 (median 
redshift z = 0.4, and 0.6 < z < 0.9 (median redshift z = 0.7). 
The two lower redshift bins comprise the majority of objects in the 
sample, as can be seen in Fig. [5] 

The centres of the bins in absolute magnitude were chosen in 
a way which ensures a fair representation of the bright end of the 
luminosity function, i.e. the bin centres at the bright end correspond 
roughly to the absolute magnitudes of the brightest objects in that 
bin. 

Fig. llOl shows the results for the if -band luminosity function 
of galaxies from spectroscopic observations in different redshift 
bins compare d to the luminosi ty f unctions determined fro m the lo- 
cal samples o flLovedavl J2OOC ) and K ochanek et all 1200 ll) .Clearly 
there is a good general agreement between the luminosity function 
measured at z = 0.2 and the ones from local galaxy samples, al- 
though our data seem to sug gest a some what smaller value for . 
Fitting a Schechter function l Schechter 1976t see also equations |5J 
and {7) for the parametrisation of the function) to the luminosity 
function yields the parameters M^- — 5 log h = (—23.79 ± 0.24) 
mag and $* K = (1.11 ± 0.12) x 10 -2 h 3 Mpc -3 . We used a 
fixed value of ape = —1.1 0, cl ose to the values deriv ed for local 
galaxies bv lLoveda^ (2000) and Koc hanek et al] 1200 lh . Contrary 
to the J-band luminosity function discussed below, our normali- 
sation also nicely agrees with the Cole et al. (2001) measurement, 
although they derive a considerably shallower faint-end slope of 
q = —0.93 ± 0.04. The errors were derived by running Monte- 
Carlo simulations with 100 000 iterations, taking into account the 
errors due to the binning in absolute magnitudes. Excluding the 



4.4 The J-band luminosity function 

The results for the rest-frame J-band luminosity function of galax- 
ies from our spectroscopic catalogue are shown in Fig. ll II Note that 
this is the first determination of the J-band luminosit y function at 
higher reds hifts following th e local measurements by Balo gh et alJ 
(200l|) and lCole et all 1200 ll) . 

The Schechter parameters (see equations J6j and for the 
definition of the Schechter function) derived in the lowest redshift 
bin (0.1 < z < 0.3, median redshift z = 0.2) are M} — 5 log h = 
(-22.45 ± 0.24) mag, $j = (1.49 ± 0.22) x 10 -2 h 3 Mpc" 3 
with the faint-end slope a.j set to a.j = —1.00. The errors were de- 
rived by running Monte-Carlo simulations with 100 000 iterations, 
taking into account the er rors due to the bin ning in absolute mag- 
nitudes. Compared to the lCole et alJ 12001) sample, our measure- 
ment seems to favour a similar characteristic magnitude M}, but 
a somewhat larger v alue of the normalisation "J?^- . Unfortunately, 
iBaloeh et alJ 120011) . the only other work on the J-band luminos- 
ity function available so far, do not derive the normalisation, thus 
the reason for the discrepancy remain s unclear. In the i f -band, the 
value of the normalisation derived bv lColeetalJl200lh is in good 
agreement with other measurements, although maybe a bit on the 
low side. 

For the intermediate redshift bin 0.3 < z < 0.6 (me- 
dian redshift z = 0.4), the Schechter parameters derived from 
our data are Mj - 5 log h = (-23.06 ± 0.24) mag, = 
(0.76 ± 0.25) h 3 Mpc -3 , where a.j was again set to a value of 
— 1.00. Thus, similar to the if -band, we see a mild evolution of the 
luminosity function with a higher characteristic luminosity and a 
smaller normalisation. The error contours for the Schechter param- 
eters of the luminosity function in the two lower redshift bins are 
shown in the same Figure. These contours were derived from the 
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Figure 9. Left panel: Model spectral energy distributions used to compute k(z) corrections (see text for details). Middle panel: i\~-band k(z) corrections as a 
function of redshift for the spectral energy distributions shown in the left panel. Right panel: Same as middle panel, but for the J band. 




Figure 10. The ii"-band luminosity function of galaxies from spectroscopic observations of MUNICS galaxies. Upper left panel: The luminosity function 
in the redshift range 0.1 < z < 0.3 (median redshift z = 0.2; filled symbols). The vertical error bars give the Poissonian error, while the binning error 
was esti mated as b/ \ZT2 with the size b of the bin. Also shown are the measurements of the local A'-band luminosity function bv lLovedavl J2000; dotted 
line) and Kochane k et alJ 1200 it dashed line), as well as a Schechter approximation to the MUNICS data (solid line). The parameters of the Schechter fit 
are M* K - 5 log h = (-23.79 ± 0.24) mag, §* K = (1.11 ± 0.12) X 10~ 2 h 3 Mpc" 3 , and a K = -1.10 (fixed). Upper right panel: The same 
as before, but in the redshift range 0.3 < z < 0.6 (median redshift z = 0.4), with Schechter parameters M* K — 51og/i = (—24.04 ± 0.26) mag, 
<S>* K = (0.71 ± 0.25) x 10~ 2 h 3 Mpc" 3 , and cik = —1.10 (fixed). Lower left panel: Error contours (la, 2<r, and 3a) for the Schechter parameters AfJ. 
and Q* K from the x 2 distribution. The filled circle indicates the average value of local measurements from Table EH the ' 'error bars" give an idea of the 
variation between different authors' measurements. Lower right panel: The luminosity function in the interval 0.6 < z < 0.9 (median redshift z = 0.70; 
filled circles), compared to the results from the two lower redshift (z = 0.2, squares; z = 0.4, triangles) as well as the local measurements. 
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Figure 11. The J-band luminosity function of galaxies from spectroscopic observations of MUNICS galaxies. Upper left panel: The luminosity function 
in the redshift range 0.1 < z < 0.3 (median redshift z = 0.2; filled symbols). The vertical error bars give the Poissonian error, w hile the binning er ror 
was estimated from b/\/l2 with the size b of the bin. Also shown are the measurement of the local J-band luminosity function bv lCole et all 1200 ll) as 
well as a Schechter approximation to the MUNICS data (solid line). The parameters of the Schechter fit are M*j — 51og/i = (—22.45 ± 0.24) mag, 
$* = (1.49 ±0.22) X 10" 2 h s Mpc -3 , and aj = -1.00 (fixed). Upper right panel: The same as before, but in the redshift range 0.3 < z < 0.6 (median 
redshift z = 0.4), with Schechter parameters M* - 5 log h = (-23.06 ± 0.24) mag, $} = (0.76 ± 0.25) h 3 Mpc -3 , and ctj = -1.00 (fixed). Lower 
left panel: Error contours (la, la, and 3a) for the Schechte r param eters Mj and $ j in the two lower redshift bins from the x 2 distribution. The filled circle 
indicates the value of the local measurement bv lCole et all l200lli with the appropriate error bar. Lower right panel: The luminosity function in the interval 
0.6 < z < 0.9 (median redshift z = 0.70; filled circles) compared to the luminosity function in the two lower redshift bins (z = 0.2, squares; z = 0.4, 
triangles) and the measurements of the local J-band luminosity function bv lCole et ali 120011) . 



X 2 distribution. At the highest redshifts probed by our sample, the 
evolution of the J-band luminosity function seems to be confirmed, 
although the statistics becomes rather poor in this case. 



5 DISCUSSION 

To quantify the redshift evolution of the near-infrared luminosity 
functions, we performed both Kolmogorov-Smirnov tests for the 
cumulative distributions of absolute magnitudes and a x 2 analysis 
for the redshift evolution of the Schechter parameters. 

5.1 Kolmogorov-Smirnov tests 

One method to compare the measurements of the near-infrared lu- 
minosity functions at various redshifts is to apply the Kolmogorov- 
Smirnov test to the cumulative distribution of objects in absolute 
magnitudes. Of course, completeness and 1/Vmax corrections have 
to be applied to each individual object entering the distribution. The 



advantage of this method is that it uses the absolute magnitudes as 
measured in the sample without binning. Furthermore, it is inde- 
pendent of the relative normalisation of the samples one wants to 
compare. 

We firstly have run this test by comparing the cumulative lu- 
minosity distributions from our data in the three redshift bins to 
the local measurements of the luminosity function. The results are 
shown in Table fTOI for the 7A"-band, and Table fTTI for the J-band, 
respectively. Furthermore, we show the resulting cumulative dis- 
tributions in Figures fl2l and [T3l for the K band and the J band, 
respectively. 

While the cumulative distributions are in fair agreement for 
the lower redshift bin (0.1 < z < 0.3), they are significantly 
different for the two higher redshift bins (0.3 < z < 0.6 and 
0.6 < z < 0.9, respectively), confirming the result found from 
the fitting of Schechter parameters M* and $*. 

As an additional test we have applied the Kolmogorov- 
Smirnov test to the cumulative distributions from the two lower 
redshift bins. The results of this test are shown in Fig. 1141 where 
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Figure 12. Resulting cumulative distribution functions for the absolute A"-band magnitudes from the MUNICS spectroscopic catalogue (solid line) as com- 
pared to the local measurements by Loveday (2000; dotted line) and Kochan e'k et alJ 1200 it dashed line). The diagram on the left-hand side compares the 
0.1 < z < 0.3 redshift bin to the local one, the middle diagram shows the same for the 0.3 < z < 0.6 redshift interval, and the diagram on the right-hand 
side compares the distributions in the 0.6 < z < 0.9 interval to the local ones. The values quoted in the diagrams give the probabilities that both cumulative 
distributions are drawn from the same population. 
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Figure 13. Resulting cumulative distribution functions for the absolute J-band magnitudes from the MUNICS spectroscopic catalogue (solid line) as compared 
to the local measurements by Cole et al. 12001; dashed line). The diagram on the left-hand side compares the 0.1 < z < 0.3 redshift bin to the local one, the 
middle diagram shows the same for the 0.3 < z < 0.6 redshift interval, and the diagram on the right-hand side compares the distribution in the 0.6 < z < 0.9 
interval to the local one. The values quoted in the plots give the probabilities that both cumulative distributions are drawn from the same population. 
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Figure 14. Cumulative distributions in absolute magnitude for the two lower redshift bins and results from a Kolmogorov-Smirnov test applied to these data, 
both for the K band (left panel) and the J band (right panel). 
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Table 8. The uncorrected values $„(Mj{) of the if -band luminosity func- 
tion, the completeness-corrected values and their errors d<E>(Mx) 
for the individual redshift intervals and magnitude bins. Here, 'complete- 
ness correction' refers to the correction of incompleteness of the photomet- 
ric catalogue, the incompleteness of the spectroscopic sample with respect 
to the photometric one, and to the V/Vm ax correction. 



Table 10. Probabilities for compatible cumulative distributions of absolute 
magnitudes derived from a Kolmogorov-Smirnov test. The comparison dis- 
tributions are derive d from the local luminosity function measurements by 
Lovedav 12000) and Kochanek et al. 1 2001). The magnitudes in the second 
column are upper limits to the range of absolute magnitudes for which the 
distribution has been computed. 



Mk — 5 log h 
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Median redshift 
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-3 
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5.69 x 10" 3 
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-3 
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7.64 x 10" 3 


9.68 x 10- 


-3 
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3 


-21.825 
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1.05 x 10" 


-2 
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-2 


4.97 x 10" 
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Median redshift 


<z> = 0.4: 
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-5 


3.90 x 10- 
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1.65 x 10" 4 


4.20 x 10- 


-4 


2.37 x 10- 
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Redshift Mk [mag] Lovedav (2000) Kochanek et al. (2001) 



Table 9. The uncorrected values $u(Mj) of the J-band luminosity func- 
tion, the completeness-corrected values <E>(Mj), and their errors d&(Mj) 
for the individual redshift intervals and magnitude bins. Here, 'complete- 
ness correction' refers to the correction of incompleteness of the photomet- 
ric catalogue, the incompleteness of the spectroscopic sample with respect 
to the photometric one, and to the V/Vm ax correction. 
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Median redshift (z) = 0.7: 
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2.24 X 10~ 4 
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(z) = 0.2 
(z) = 0.4 
(z) = 0.7 



-21.5 
-23.0 
-24.0 



35.73 % 
0.88 % 
1.61 % 



26.10% 
0.15 % 
0.37 % 



Table 11. Probabilities for compatible cumulative distributions of abso- 
lute magnitudes derived from a Kolmogorov-Smirnov test. The comparison 
distribution is derived from the local luminosity function measurement by 
IColeetaljEool . The magnitudes in the second column are upper limits to 
the range of absolute magnitudes for which the distribution has been com- 
puted. 



Redshift 


Mj [mag] 


Cole et al. C2001) 


(z) = 0.2 


-20.3 


21.51 % 


(z) = 0.4 


-21.6 


0.01 % 


(z) = 0.7 


-22.9 


0.07 % 



one can see that the probabilities that both distributions are drawn 
from the same parent population are very small. 



5.2 Likelihood analysis of luminosity function evolution 

In this Section we will describe a \ 2 analysis of the redshift evo- 
lution of the near-infrared luminosity functions. This test uses the 
Schechter parametrisation 



(6) 



of the luminosity function, where L* is the characteristic lu- 
minosity, a the fain t-end slope, and <E>* the normalisation of the 
luminosity function <Schechterll97r3) . The corresponding equation 
in absolute magnitudes reads 



* (M) = - $* In 10 10' 

5 



exp 



-10 



0.4(M* -M)(l+a) 



(7) 



To estimate the rate of evolution of the parameters with red- 
shift, we define evolution parameters /i and v as follows: 



(z) = (0) (1 + /w) , 
M* (z) = M* (0) + vz, and 
a (z) — a (0) = a. 



(8) 



Note that the faint end of the luminosity function cannot be 
determined very well from our data, thus we leave the faint-end 
slope a of the Schechter luminosity function fixed, as we have also 
done during the fitting of a Schechter function to our data. 
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To quantify the redshift evolution of $* and M* we now com- 
pare our luminosity function data in all redshift bins with the local 
Schechter function evolved according to equation |8j to the appro- 
priate redshift. We do this for a grid of values for /i and v, and 
calculate the value of \ 2 for eacn grid point according to 

N o 

2, lv MM^^-^jM^^v, Zi )Y 

x ^' V > ~ n Z^i a 2 ' 

i — 1 1 

where <!>(M, z) is the measurement of the luminosity func- 
tion at median redshift z in the magnitude bin centred on M, 
ty(M, /x, u, z) is the local Schechter function evolved according to 
the evolution model defined in equation J8} to the redshift z, en is 
the RMS error of the luminosity function value, and n is the num- 
ber of free parameters of the approximation, i.e. the number of data 
points used minus the number of parameters derived from the fit- 
ting. 

We want to compare our measurement of the A-band lu- 
minosity function with the Schechter approximatio ns to the lo- 
cal determinations. We use the measurements by lLoved"avl l2000l) 
and K ochanek et alJ 1200 lh for the if -band (since the luminosity 
function parameters derived from local samples are very similar 
anyway), a nd the local J-band luminosity function is the one by 
ICole et ail j200lh. The Schechter parameters derived by those au- 
thors are shown in Tables IT2l and 1 1 3l Chosing a shall ower faint- 
end sl ope in the A band, similar to the one derived bv lCole et alJ 
1200 ll) . changes the result slightly, but - within the errors - not 
significantly. 

To avoid that data points with large completeness correction 
factors affect the result, we exclude all luminosity function mea- 
surements with a total correction factor (photometric incomplete- 
ness, spectroscopic incompleteness, and U/V max correction) larger 
than three. 

The result of the likelihood analysis is shown in Fig. 1151 
For the A"-band, we compare our measurements at redshifts 0.2, 
0.4, and 0.7 to the lo cal measurements by iLovedavl l200fj) and 
iKochanek et alJ 1200 lh . and to the average of their Schechter pa- 
rameters. We detect a brightening of AM^/Az ~ 0.70 magni- 
tudes, and a decline of the number density of objects to redshift 
one. The decrease of <£>| c with redshift is obviously quite strongly 
dependent on the parameters of the local luminosity function, how- 
ever, for the average value we derive A$* K / [Q* K Az) ~ —0.35. 
These results also give quantitative estimates of the evolution which 
can already be seen in the Schechter p arameters der i ved fr om our 
data, see Table|8|for details. Note that lHuang et alJ i2002^ derive 
a significantly brighter M* , a slightly larger normalisation, and a 
steeper faint-end slope, which they ascribe to redshift selection ef- 
fects. If their measurements are valid, the brightening to redshift 
one would be smaller, whereas the evolution in number density 
would be even larger. 

Within the errors, the results found in this work agree nicely 
with the measurements of the A-band luminosity function derived 
from the full MUNICS sample based on photometric redshifts. First 
results are shown in MUNICS III and show the same trend for the 
evolution of the luminosity function with redshift. A more detailed 
analysis will be shown in MUNICS II, where the evolutionary trend 
with a brightening of 0.5 to 0.7 mag and a decrease in number den- 
sity of roughly 25 per cent to redshift one is confirmed. 

In the case of the J-band (right panel of Fig. 1151 . the evolu- 
tion of the luminosity function is obviously not constrained. This 
is also apparent from the error contours of the Schechter parame- 
ters shown in Fig. lllK lower left panel), where one can see that the 



local measurement bv lCole et alJ 1200 ll) has a characteristic mag- 
nitude similar to the one derived in our lowest redshift bin, but a 
normalisation in between the ones derived from our two lower red- 
shift intervals, thus making any conclusions about evolution with 
respect to the local sample difficult. Nevertheless, we note that our 
J-band luminosity function data seem to confirm the trend seen for 
the A band, which is also evident from the Kolmogorov-Smimov 
tests presented in Section l5~T1 



6 CONCLUSIONS 

We have presented spectroscopic follow-up observations of galax- 
ies from the Munich Near-Infrared Cluster Survey (MUNICS), de- 
scribed the observations, the data-reduction and the properties of 
the spectroscopic sample. Furthermore we have presented the rest- 
frame A-band luminosity function for galaxies at median redshifts 
of z — 0.2, z — 0.4, and z = 0.7. The Schechter parameters de- 
rived at redshift z = 0.2 are M£-5 log h = (-23.79±0.24) mag, 
<S>K = (1.11 ± 0.12) x 10" 2 h 3 Mpc -3 for fixed a K = -1.10 
as measured locally, in good agreement to values derived at low 
redshifts. At redshift z = 0.4, however, the Schechter parame- 
ters are M* K - 5 log h = (-24.04 ± 0.26) mag and ®* K = 
(0.71 ±0.25) x 10" 2 h 3 Mpc~ 3 (for the same value of a K ). Thus 
the value for the characteristic luminosity is somewhat larger and 
the normalisation smaller than the ones derived locally. This is con- 
firmed by Kolmogorov-Smirnov tests for the cumulative distribu- 
tions in absolute magnitudes, and by a x 2 analysis for the evolution 
of the luminosity function. From the latter we find mild evolution 
in magnitudes (AMj { = —0.70 ± 0.30 mag) and number densi- 
ties (A$* K /<b* K = -0.35 ± 0.30) to redshift one. Furthermore, 
we have presented the first measurement of the J-band luminosity 
function of galaxies at higher redshifts with Schechter parameters 
M} - 5 log h = (-22.45 ± 0.24) mag, $} = (1.49 ± 0.22) x 
10~ 2 /i 3 Mpc" 3 forz = 0.2, and M} -5 log /t = (-23.06±0.24) 
mag, $} = (0.76 ± 0.25) h? Mpc" 3 for z = 0.4, showing the 
same trend of evolution as the A-band luminosity function. The 
faint-end slope a j set to a value of — 1 .00 in both redshift bins. The 
evolutionary trend of the near-infrared luminosity of A-selected 
galaxies described above is consistent with expectations from pure 
luminosity evolution, while the decrease in number density with 
redshift can be understood in the context of hierarchical galaxy for- 
mation models. 
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Figure 15. Left panel: Result of the estimation of the redshift evolution of th e Schechter parameters and <fr* K between z = 0.4 (MUNICS luminosity 
function) and the local universe from lLovedavl <200(j dotted line), iKochanek et"aH feOOlt dashed line), and for an average of these two local measurements 
(solid line) as derived from a x 2 approach (see text for details). The contours correspond to la, 2<r, and 3c confidence level. The dQ* /dz = and 
dM* /dz = li nes in dicate the non-evolution values. Right panel: The same for the J-band luminosity function. In this case, the local measurement is taken 
from Cole et al. 1 2001, dotted line), with the appropriate Schechter parameters from Table ll3l and the evolution is very badly constrained. 



Table 12. Summary of samples used for deriving the i-sT-band field-galaxy luminosity function, giving the source for the data, the approximate limiting 
magnitude (either in the optical or in the near-infrared), the approximate area, the number of objects, the redshift range of the survey (either as interval or 
as median redshift) and the parameters of the Schechter function derived from the data. Parameters without error estimate were kept fixed during the fitting 
proce dure. Where necessary, the values quoted have been converted from the cosmology originally used to the f2 m = 0.3, Ha = 0.7 cosmology following 
IColeetaljkOQll) . 



Source 


Limit 
[mag] 


Area 
[arcmin 2 ] 


# 


2 


M* K - 5 log h 
[mag] 


a K 


$* 

if 


Mobasher et al. (T9931 1 


Bj < 17, K <13.0 




181 


0.0... 0.1 


-23.37 ± 0.30 


-1.00 ± 0.30 


1.12 ±0.16 


Songaila et al. (T9941 2 


K < 14.5... 20.0 


5544. . . 5 












Glazebrook et al. (1995) 3 


K <17.0 


552 


124 


0.0... 0.8 


-23.14 ± 0.23 


-1.04 


2.22 ± 0.53 


Cowie et al. (19961 4 


K <19.5 


26 


254 


0.2... 1.0 


-23.49 ±0.10 


-1.25 ± 0.15 


0.80 ± 0.20 


Gardner et al. ( 19971 5 


K <15.0 


15840 


567 


0.14 


-23.30 ± 0.17 


-1.00 ±0.24 


1.44 ±0.20 


Szokolvetal. a998) 6 


R F < 18.5, K <16.5 


2160 


175 


0.0... 0.4 


-23.80 ± 0.30 


-1.30 ±0.20 


0.86 ± 0.29 


Lovedav (20001 7 


K <12.0 




363 


0.05 


-23.58 ± 0.42 


— 1.16 ± 0.19 


1.20 ±0.08 


Kochanek et al. (20011 8 


K 20 < 11.25 




4192 


0.02 


-23.43 ± 0.05 


-1.09 ±0.06 


1.16 ±0.10 


Cole et al. (20011 9 


bj < 19.5, K < 13 


2.2 10 6 


17173 


0.05 


-23.36 ± 0.02 


-0.93 ± 0.04 


1.16 ±0.17 


Baloghetal. (2001 I 10 


Be < 15 . . . 17, K < 13 






0.0... 0.18 


-23.48 ± 0.08 


— 1.10 ± 0.14 




Huang et al. f2002) n 


K < 15 




1065 


0.138 


-23.70 ± 0.08 


-1.39 ±0.09 


1.30 ±0.20 


MUNICS (this work) 


K < 17.5... 19.0 


649 


157 
145 


0.1... 0.3 
0.3... 0.6 


-23.79 ± 0.24 
-24.04 ± 0.26 


-1.10 
-1.10 


1.11 ± 0.12 
0.71 ± 0.25 



1 Based on the optically-selected Anglo-Australian Redshift S urvey {P eterson et al|l986ll. A correction of +0.22 magnitudes is sometimes applied to their 
result because of their method of calculating k(z) corrections iGlazebrook et alll995t) 

2 The area of the Songail a et all 1 1 99411 sample actually continuously decreases with K magnitude. The values in the table are intende d to give a rough 
impression of the parameters of the catalogue. No luminosity function is derived by the authors, however, the sample has been used bv lCowie* etal. (1996) in 
combination with their deep sample to derive the i-sT-band luminosity function in various redshift intervals. 

3 A correction of —0 .30 magnitudes is often applied to their result due to the use of fixed-aperture photometry i Glazebro ok et alll995t) . 
JJCoujcjjUiJ^Wf: I determine the luminosity function in four redshift bins from a combination of their deep sample with the shallower samples from 
Songail a et al J 1 1 9941) . The values for the Schechter parameters given in the table are from their fit to the luminosity function over the whole redshift range, 
with errors estimated from the dispersion of the values for different redshift intervals. 

5 The spectroscopic catalogue o f lGardneret ail 1 19971) is sparsely selected on the large photometric sample with the parameters given in the table due to 
geometric limitations imposed by the fibre spectrograph used for the observations. ^^^^^^^^^^ 

B Based on redshift from the optically-selected Kitt-Peak Galaxy Redshift Survey <Munn et alll997t). 

Based on ii'-band imaging of galaxies from the optically-selected Stromlo-APM galaxy survey i Mad dox et alll99CHLovedav et all 1996). 
x Based on a combination of the 2MASS catalogue |jm£tteUij2000l) with the CfA2 i Geller & Huchr jl9891) and UZC <Falco et alll999l) surveys. 

9 Based on a combination of the 2MASS catalogue tlangt^e^JgggfA with the 2dFGRS fFolkes et alll999l) . 

10 Based on a combination of the 2MASS catalogue Llarrett et all20 ()0 ) with the Las Campanas Redshift Survey i Shectma n et alll996l) . 

11 Hawaii-AAO K-band Galaxy Redshift Survey. 
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Table 13. Summary of samples used for deriving the J-band field-galaxy luminosity function, giving the source for the data, the approximate limiting 
magnitude (either in the optical or in the near-infrared), the approximate area, the number of objects, the redshift range of the survey (either as interval or 
as median redshift), and the parameters of the Schechter function derived from the data. Parameters without error estimate were kept fixed during the fitting 
procedure. 



Source 


Limit 
[mag] 


Area 
[arcmin 2 ] 


# 


2 


M* — 5 log h 
[mag] 


aj 


[10- 2 /i 3 Mpc- 3 ] 


Cole et al. (2001) 1 
Balogh et al. (2001) 2 


bj < 19.5, K < 13 
Rc < 15 . . . 17, K < 13 


2.2 10 6 


17173 


0.05 
0.0... 0.18 


-22.36 ±0.02 
-22.23 ± 0.07 


-0.93 ± 0.04 
-0.96 ±0.12 


1.08 ±0.16 


MUNICS (this work) 


K < 17.5. . . 19.0 


649 


132 
145 


0.1... 0.3 
0.3... 0.6 


-22.45 ± 0.24 
-23.06 ± 0.24 


-1.00 
-1.00 


1.49 ± 0.22 
0.76 ± 0.25 



1 Based on a combination of the 2MASS catalogue 

2 Based on a combination of the 2MASS catalogue 



Jarrett et al|200( I with the 2dFGRS iFolkes et alll999h . 
Jarrett et all20 00) with the Las Campanas Redshift Survey IShectman et alll99d 
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